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We report on our study of light hadron spectrum and quark masses in QCD with two flavors of dynamical 
quarks. Simulations are made with the plaquette gauge action and the non-perturbatively O(a) improved Wilson 
quark action. We simulate 5 sea qaurk masses corresponding to mps / mv — 0.8-0.6 at /3 = 5.2 on 12 3 x48, 16 3 x48 
and 20 s x 48 lattices. A comparison with previous calculations in quenched QCD indicates sea quark effects in 
meson and quark masses. 



1. Introduction 

One of the major goals of lattice QCD simu- 
lation is to confirm the validity of QCD as the 
theory of strong interaction in the low energy re- 
gion by comparing its prediction for the hadron 
spectrum to experiment. Because of huge compu- 
tational demand to perform simulations with dy- 
namical quarks, however, many works have been 
forced to neglect effects of dynamical sea quarks. 
In this quenched approximation, the CP-PACS 
collaboration found some deviation in the light 
hadron spectrum from experiment [Q, which mo- 
tivates us to investigate sea quark effects by per- 
forming simulations of QCD with the realistic 
number of dynamical quark flavors. 

The JLQCD collaboration started numerical 
simulations of dynamical QCD on a supercom- 
puter Hitachi SR8000 model Fl, which is newly 
installed at KEK on March 2000. It has a 1.2 
TFlops of peak performance and provides about 
400 sustained GFlops for our simulation code. We 
carry out simulations of QCD with two flavors 
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of dynamical quarks and investigate the chiral 
extrapolation and finite size effects in the light 
hadron spectrum, as a step toward studies of 
QCD with 2+1 flavors (i.e. u, d plus s quarks). 

2. Simulations 

We study QCD with two degenerate flavors of 
dynamical quarks, which are identified with u 
and d quarks; the strange quark is treated in the 
quenched approximation. We employ the O(a) 
improved quark action^] with the clover coeffi- 
cient csw determined non-perturbatively by the 
ALPHA collaboration^. Simulations are per- 
formed at (3 = 5.2 and csw = 2.02, for which 
a -1 evaluated at the physical mass of dynami- 
cal ud quarks is approximately 2 GeV and hence 
scaling violation is not expected to be too large. 
We simulate 5 sea quark masses in the range 
mps/mv — 0.8-0.6. Since finite size effects could 
be more important with dynamical quarks, we 
perform simulations on three lattices with differ- 
ent spatial sizes, 12 3 x 48, 16 3 x 48 and 20 3 x 48. 
We have accumulated 3000 thermalizcd HMC tra- 



Table 1 



Simulation parameters. We also list the lattice 
spacing a ra fixed by r = 0.49 fm for 16 3 x 48 
lattices. 



lattice 




#traj. 


mps/mv 


a r „ [fm] 


16 3 x48 


0.1340 


3000 


0.802(4) 


0.1288(10) 




0.1343 


3000 


0.781(6) 


0.1201(10) 




0.1346 


3000 


0.743(6) 


0.1127(7) 




0.1350 


3000 


0.714(9) 


0.1084(7) 




0.1355 


3000 


0.596(18) 


0.1006(5) 


12 3 x48 


0.1346 


3000 


0.735(9) 






0.1350 


3000 


0.695(15) 




20 3 x48 


0.1346 


2000 


0.756(8) 






0.1350 


2800 


0.706(9) 





jectories on the 16 3 x 48 lattice, while simulations 
on other lattices are in progress. Other details of 
simulation parameters are summarized in Table [|. 

3. Static quark potential 

We calculate the static quark potential on the 
16 3 x 48 lattice using the smeared Wilson loops. 
Since we do not observe any signal of flattening 
of the potential due to string breaking, potential 
data are parameterized with the form 

V(v) = V -a/\r\-g-SV(r)+a-\r\, (1) 

where 5V(r) = (G(r) — l/|r|) represents a correc- 
tion to the short-distance Coulomb term calcu- 
lated with the lattice gluon propagator^]. Fig- 
ure [I] shows a plot obtained for the heaviest sea 
quark mass, where we find that deviation of the 
potential data from the fit curve is at most 1% 
in the fitting range of r. We determine the Som- 
mer scale tq j(| from the parameterization (Q). 
The lattice spacing determined with the condi- 
tion r = 0.49 fm is listed in Table |l|. 

4. Light meson mass measurement 

We calculate light hadron correlators for 6 val- 
ues of the valence quark masses in the range 
mps/my ~ 0.80-0.50 for each sea quark mass. 
Measurements have been completed for the 16 3 x 
48 lattice, while only the degenerate case {K va \ = 
K sca ) is measured on 12 3 x 48 and 20 3 x 48. Most 
of the analysis results discussed in this talk are 
based on the 16 3 x 48 lattice data except for the 
finite size effect in the following. 




r/a 

Figure 1 . Corrected static quark potential V(r)+ 
g ■ SV(r) at if sca = 0.1340 on 16 3 x 48 lattice. The 
solid line represents the fitting curve (El). 



Figure |2j shows the effective mass of degenerate 
(^vai = -Ksea) pseudo-scalar meson at the second 
lightest sea quark mass. Data are shown for three 
different spatial volumes 12 3 , 16 3 and 20 3 . We 
find that the fitted masses on the two larger lat- 
tices are in good agreement with each other both 
for the pseudo-scalar and vector mesons, as shown 
in Figure |^. This is also observed for other lattices 
with heavier sea quarks. This suggests that finite 
size effects are already small for the 16 3 lattice. 
We are currently extending similar measurements 
to the lightest sea quark K sca = 0.1355, for which 
finite size effects are expected to be more impor- 
tant. 

5. Sea quark effect in meson masses 

Figure |] shows the vector meson mass as 
a function of the pseudo-scalar meson mass 
squared. Since the effective lattice spacing de- 
creases as the sea quark mass becomes smaller, 
which complicates the chiral extrapolation if 
made using the lattice unit, we plot the meson 
masses normalized by a physical quantity ro mea- 
sured for each sea quark mass. We also plot the 
quenched data @ in Figure |J, where we clearly 
see an indication of sea quark effect, i.e. slope in 
two-flavor QCD is significantly larger than that 
in quenched QCD. This leads to a larger hyper- 
fine splitting of the strange meson as we shall see 
later. We also find that our results at different sea 
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Figure 2. Effective masses for the pseudo-scalar 
meson at K sca _ — 0.f350 on three different vol- 
umes. Solid lines represent the central value of 
the fit result. 

quark masses lie almost on one curve. This indi- 
cates that the sea quark mass dependence in this 
quantity is not large in our simulated region of 
-Kgoa and more precise calculation is needed, par- 
ticularly for the vector mesons, to see it clearly. 

6. Chiral extrapolation 

We employ two different strategies for the chi- 
ral extrapolation. In our main strategy, which we 
call method (A), we use meson and quark masses 
normalized by ro throughout the analysis. For 
instance, the pseudo-scalar meson mass squared 
is fit to the form 



(m PS ro) 2 = A s (m sea ,r ) + A v {m va \r ), 



(2) 



where m aca , and m va i are respectively sea and av- 
erage valence quark masses defined through the 
vector Ward identity (VWI) relation to 9i vwi = 
(1/K - l/K c )/2. With this method, data is de- 
scribed very well by a simple linear fit. In the 
other strategy, method (B), we make the chiral 
extrapolation using meson and quark masses in 
the lattice unit. In this case, significant curvature 
is found in fits for both pseudo-scalar and vector 
meson masses, and therefore we have to introduce 
quadratic terms in the fitting function. We find 
that both fits give consistent results, although the 
fit error is generally larger for the method (B). 

In each method, the physical quark masses m u d 
and m s are determined by tuning meson mass 
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Figure 3. Fit results of pseudo-scalar (circles) 
and vector (squares) meson masses as a function 
of (16/L) 3 . Results on V — I6 3 are shown also 
with dashed lines for a guide of eye. 



ratios to their experimental values. We use m„ 
and trip to fix m uc i and the lattice spacing. The 
strange quark mass m s is determined from either 
m^/wp or m^/rrip. 

Our results for the strange vector meson masses 
(K* and </>), with pseudo-scalar (K) used to fix 
the strange quark mass, are shown in Figure |^ 
by filled circles (method (A)) or by filled squares 
(method (B)). For comparison, we also plot the 
CP-PACS data in the quenched approximation (!]] 
and the recent two-flavor result || . A trend that 
the sea quark effect pushes up the vector meson 
masses toward the experimental values is seen in 
this plot. 



7. Strange quark mass 

The bare strange quark mass is calculated us- 
ing either of two definitions, i.e. one from VWI 
and the other from the axial vector Ward identity 
(AWI), 2m,, AW i = (a^4(*)P(0)t)/(P(t)P(0)t>, 
where A^t) and P(t) are the axial vector current 
and pseudo-scalar density. We use the improved 



axial vector current A 



(imp) 



(t)=A 



l {t)+c A d ti P{t), 
with ca calculated at one-loop g. The contin- 
uum MS quark mass is obtained using one-loop 
matching |||| at scale fi = 1/a and evolved to 
H = 2 GeV with 3-loop (3 function 0. 

Our results for the strange quark mass using 
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Figure 4. (mv r o) vs (tops^o) 2 m two flavor 
QCD (filled symbols). Open symbols represent 
results in quenched QCD in Refs. @. Experimen- 
tal values are plotted with open diamonds using 
r = 0.49 fin. 

the method (A) are 

94(2) MeV (VWI), 



m s (2GeV) - ^ gg ^ Mey (AW ^ 
with K used as input, or 



m s (2GeV) = 



109(4) MeV (VWI), 
102(6) MeV (AWI), 



(3) 



(4) 



with <f> used as input. The method (B) gives con- 
sistent results. 

Only the statistical errors are shown in (g) and 
(^). The systematic error is more significant, as 
indicated by the disagreement between VWI and 
AWI, and the difference between (|J) and (||), al- 
though the latter is smaller than in the quenched 
results jlj. One of the most important sources of 
the systematic error is the use of the one-loop 
pcrturbative Z factor, for which a naive order 
counting gives 0(5%). In addition, the discretiza- 
tion error of 0(a 2 ) or the quenching effect of the 
strange sea quark could also be important. A re- 
cent two-flavor simulation of CP-PACS predicted 
m s (2GeV) = 88tg MeV or 90±^ MeV in the con- 
tinuum limit with the K or <f> used as input ||. 
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Figure 5. Strange meson masses with K- input. 
Filled circles and squares are our full QCD results 
using method (A) and (B), respectively. Open 
triangles and squares are CP-PACS's result in 
quenched]!]] and full QCDQ]. Horizontal lines 
show experimental values. 
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